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We investigate the bulk band structures and the surface states of Bi2Se3 and Bi2Te3 topological 
insulators using first-principles many-body perturbation theory based on the GW approximation. 
The quasiparticle self-energy corrections introduce significant changes to the bulk band structures, 
while their effect on the band gaps is opposite in the band-inversion regime compared to the usual 
situation without band inversion. Parametrized "scissors operators" derived from the bulk studies 
are then used to investigate the electronic structure of slab models which exhibit topologically 
protected surface states. The results including self-energy corrections reveal significant shifts of 
the Dirac points relative to the bulk bands and large gap openings resulting from the interactions 
between the surface states across the thin slab, both in agreement with experimental data. 
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The recently discovered bulk topological insulators 
(TIs), a new class of semiconducting materials charac- 
terized by the presence of spin-helical surface states re- 
sulting from strong spin-orbit (SO) interactions, have 
quickly become a subject of intense research It is 
believed that TIs are actually not uncommon among the 
heavy-element materials. Owing to their novel electronic 
properties, such as the suppression of backscattering and 
intrinsic spin polarization of the surface-state charge car- 
riers, TIs are expected to find applications in the future, 
including information technology, spintronics and quan- 
tum computing. 

First-principles electronic structure calculations can 
play an important role in exploring the properties of 
known TIs as well as in guiding the search of novel 
materials. The widely used density functional theory 
(DFT) within the Kohn-Sham formalism, a workhorse 
first-principles method of condensed matter physics, has 
proved its value already in the discovery of the "second- 
generation" bulk TIs, Bi 2 Se 3 and Bi 2 Te 3 HQ. A large 
number of TIs has been predicted using this technique 
and some of these predictions have been confirmed 
experimentally However, it is broadly recognized 

that the Kohn-Sham eigenvalues of DFT fail to describe 
accurately quasiparticle energies and band gaps the 
critical properties of TIs. A recent work has highlighted 
the limitations of standard DFT approach in describing 
the topological nature of several borderline compounds 
Q. Many-body perturbation theory techniques, such as 
the GW approximation, greatly improve the accuracy of 
predicting these excited-state properties d, 0] . 

In this Letter, we investigate the self-energy effects in 
the quasiparticle bulk band structures and the surface- 
state dispersion of the reference TIs, Bi2Se3 and Bi2Te3, 
by using the first-principles GW method @. We find 
that the effects of quasiparticle corrections on the band 



structures are substantial, and show novel features result- 
ing from the interplay with SO interactions in the band- 
inversion regime. By introducing parametrized energy- 
dependent "scissors operators" based on the bulk cal- 
culations, we have extended our study to the electronic 
structure of slab models. Our application of the pro- 
posed technique to slab models which exhibit topologi- 
cally protected surface states finds significant shifts of the 
Dirac point energies relative to the bulk bands and larger 
surface-state gap openings resulting from the interactions 
between surface states across the slab. The proposed sim- 
ple corrections basically eliminate the DFT eigenvalue 
problems and yield agreement with experimental data. 

The DFT calculations were performed within the local 
density approximation (LDA) employing the Quantum- 
ESPRESSO package [lOj. We used norm-conserving 



pseudopotentials [llj and a plane-wave kinetic energy 
cutoff of 35 Ry for the wavefunctions. The quasiparticle 
energies were evaluated within the GqWq approximation 
to the electron self-energy starting from LDA results as a 
mean-field solution using the approach of Hybertsen and 
Louie [§]. The static dielectric function was calculated 
using a 10 Ry plane- wave cutoff, unoccupied bands up to 
5 Ry above the Fermi level, and extended to finite fre- 
quencies with the generalized plasmon-pole model. This 
first-principles GW methodology is implemented in the 



BerkeleyGW code [12(. Spin-orbit interactions were 



included on the final stage using the SO Hamiltonian ma- 
trix Hso(k) evaluated employing an approach described 
in Ref. [l3| in the basis of eigenfunctions of "Ho(k). That 
is, in the full two-component Hamiltonian 



«(k) = ^o(k) + ^ so (k) 



(1) 



7^0 (k) is a diagonal matrix with matrix elements being 
either LDA or GW eigenvalues in the absence of SO in- 
teractions, while %so(k) introduces off-diagonal matrix 
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FIG. 1: (color online) . Quasiparticle self-energy corrections as 
a function of LDA energies for bulk (a) Bi2Se3 and (b) Bi2Te3 
calculated without taking into account SO interactions. The 
positions of VBM and CBM are indicated by the larger filled 
and open circles, respectively. The lines correspond to the 
fitted "scissors operators". (c,d) Same plots after the SO 
matrix elements were taken into account. Note the changes 
in GW corrections calculated for the states which correspond 
to VBM and CBM in (a) and (b). 



elements. Notably, we find that the results obtained us- 
ing this method are in agreement with an explicit ap- 
proach in which GW calculations are performed starting 
from two-component LDA wavefunctions after SO inter- 
actions were taken into account (and thus, band inversion 
is present at the T point) [13]. We used experimental lat- 
tice parameters for both the bulk crystal as well as (111) 
slab models of different thickness. 

We start our discussion by considering the GW quasi- 
particle energy corrections A£ , Qp(n,k) = Ecw{n,\i) — 
Elda{h, k) for bulk Bi 2 Se 3 and Bi 2 Te 3 . Figs. QTa,b) 
show AEqp(ti, k) evaluated on a 6x6x6 fc-point grid as 
a function of LDA energy with no SO interactions taken 
into account. For convenience, we set the energies cor- 
responding to the valence band maximum (VBM) as a 
reference (i.e. A£ QP (VBM) = £ LD a(VBM) = for the 
case with no SO interactions included). When SO inter- 
actions are neglected both materials are direct band gap 
semiconductors with VBM and conduction band min- 
imum (CBM) located at the T point [dotted lines in 
Figs. Hfa,b)]. The inclusion of GW self-energy correc- 
tions increases the LDA direct gaps at T point (which 
are 0.151 eV and 0.188 eV for Bi 2 Se 3 and Bi 2 Te 3 , re- 
spectively) by 0.212 eV and 0.277 eV, respectively. Af- 
ter SO interactions have been introduced in both LDA 
and GW calculations, the values of AEqp(ti, k) barely 
change except for those which correspond to VBM and 
CBM in Figs.[TJa,b). These states, shown as large circles 



in Figure [T] change their order in energy. 

The observed seemingly counterintuitive behavior is a 
direct consequence of band inversion due to SO interac- 
tions. It can be illustrated using the k ■ p Hamiltonian 
for Bi 2 Se 3 -type materials [3j: 



n k . p (k) = 



( A4(k) 

A z k z 
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(2) 

with k± = k x ± iky and k = |k|. Without loss of gen- 
erality we assume A4(k) = — A g /2 — k 2 /2m* with a sin- 
gle m* parametrizing both the valence and conduction 
bands, eo(k) = 0, and A xy = A z = A. As the param- 
eter A g decreases, a band inversion takes place around 
k = for A g < and the band gap closes in the ab- 
sence of off-diagonal matrix elements (dotted lines in 
Fig. [3]). However, these off-diagonal SO matrix elements 
ensure non-zero band gap even in the band-inversion 
regime. For a small magnitude A g < the bands re- 
main parabolic and the gap is — A g . Further decrease of 
A g leads to the "camelback" shaped bands and the band 
gap is Aso 2Ak. The "camelback" feature is clearly 
seen for the valence bands in the LDA band structures 
of both Bi 2 Se 3 and Bi 2 Te 3 [dashed lines in Fig. [U[a,b)]. 
The physical effect incorporated in the quasiparticle self- 
energy correction is the increase of A g which is typically 
underestimated in DFT. Upon an increase of the value 
A g in the A g < regime, the energy of the valence band 
at k = increases while the energy of the conduction 
band at k = decreases. This behavior is the oppo- 
site to the "normal" situation where no band inversion 
takes place. In other words, the quasiparticle self-energy 




FIG. 2: (color online). Band structures of bulk (a) Bi2Se3 
and (b) Bi2Te3 calculated using the following theories: LDA 
with no SO (dotted lines), LDA with SO effects taken into 
account (dashed lines) and GW with SO interactions (solid 
lines). (c,d) Degree of band inversion for the valence band 
(Eq. [3]) of Bi2Se3 and Bi2Te3, respectively, calculated along 
the same fc-path as in panels (a) and (b). 
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FIG. 3: (color online). Evolution of the band structure (solid 
lines) calculated for model Hamiltonian given by Eq. ((2]) upon 
the decrease of parameter A g (see text). Dotted lines repre- 
sent the solution calculated in the absence of off-diagonal (SO) 
matrix elements. 



corrections to the inverted bands reduce the direct Kohn- 
Sham DFT gap at T. This is exactly what is observed in 
Figs. QIc,d) (large symbols) and in Figs.[2ja,b) at T. 

The band structures of bulk Bi2Se3 and Bi2Te3 includ- 
ing both SO interactions and the GW quasiparticle self- 
energy corrections [solid lines in Fig. H{a,b)] have been 
calculated using the Wannier interpolations technique 
[IE>[l(|. Interestingly, the LDA band gap of bulk Bi 2 Sc 3 
(0.29 eV) barely changes after the inclusion of GW cor- 
rections (0.30 eV). However, its character changes from 
indirect to direct in agreement with recent experiments 
17| . The surprising accuracy of LDA in predicting 
the magnitude of minimum band gap is fortuitous. In 
Bi 2 Te3, the LDA and GW indirect band gaps are 0.09 eV 
and 0.17 eV [l8j], with the latter being in good agreement 
with experiment One noticeable effect of the GW 

corrections on band dispersion is a considerable diminu- 
tion of the dip in the valence bands at the V point. This 
behavior is also consistent with the discussed two-band 
model [Eq. ©; Fig. [3]. 

In order to gain further understanding of the effects of 
GW corrections we define the degree of band inversion 
for the valence band 



BI VB (k)= J2 4( k K(k), (3) 

i— VB;j£Eunocc. 

where the eigenfunctions ipi(k) — J2j a i?(k)^j (k) °f 
Hamiltonian H(k) = "Ho(k) + 'Hso(k) are expressed in 
terms of ^(k), the eigenfunctions of 'Ho(k) which does 
not include SO interactions. The results for both LDA 
and GW methods are plotted in Figs. [2Jc,d). For both 
Bi2Se3 and Bi2Te3 band inversion takes place only in a 
limited region of the Brillouin zone around the T point 
where BIvb achieves almost 100%. The introduction of 
GW shifts somewhat reduces the extension of this region 
of band inversion, but at the T point it is still complete. 
A very similar picture was obtained for the degree of band 
inversion of the conduction band (not shown) . 

The crux of our study is, of course, to investigate the 
effects of GW quasiparticle self-energy corrections on the 



topological surface states. Addressing this problem in a 
straightforward way would require performing GW cal- 
culations for two-dimensional slab models. At present, it 
is computationally too demanding to perform converged 
GW calculations on systems of this size. To overcome 
this difficulty, we parametrize the quasiparticle correc- 
tions AEqp in terms of energy- dependent "scissors op- 
erators" AEqp = ai?LDA + bE^DA + c using the results 
of our GW calculations for bulk materials with no SO 
interactions included. Valence and conduction bands are 
fitted separately. Additionally, we require our "scissors 
operators" to reproduce exactly the quasiparticle shifts of 
the VBM and CBM, the critical components in our con- 
sideration. The fitted functions are shown in Figs.QJajb), 
and the corresponding parameters are given in Ref. l2(t 

Figures |4|a,b) show the band structures in the vicinity 
of the r point computed for (111) slabs of 5-quintuplc- 
layers (5QL) thickness using plain LDA and after ap- 
plying the proposed generalized "scissors operator" tech- 
nique. For the latter, we apply the "scissors operators" 
to the LDA Hamiltonian without SO interactions in the 
Bloch-state basis to obtain the quasiparticle Hamilto- 
nian. The SO coupling terms are then added to it and 
the Hamiltonian is diagonalized to obtain the final quasi- 
particle energies. Both methods give rise to topologi- 
cal^ protected surface states appearing as characteristic 
"Dirac cone" features at the T point but with a gap Eg 
owing to the hybridization of the surfaces states at the 



opposite surfaces of thin slabs [2l|, |22j. The inclusion 



of quasiparticle corrections results in the following two 
important changes. First, the magnitudes of E^ are en- 
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FIG. 4: (color online). Band structures calculated for 5QL 
slabs of (a) Bi 2 Se 3 and (b) Bi 2 Te 3 using LDA and GW (SO 
interactions included in both cases). (c,d) Band gaps at the V 
point as a function of slab thickness for Bi2Se3 and Bi2Te3, re- 
spectively. Experimental results are reproduced from Ref. l23l . 
The zero of energy is set at 15lda(VBM) with no SO. 
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larged, especially in Bi 2 Se 3 [Figs. 0|c,d)]. While LDA 
predicts essentially zero values of Eg for slabs thicker 
than 3QL, the dependence turns into a slow 1 /width de- 
cay after the "scissors operator" was applied [Fig. Be)]. 
This behavior is actually consistent with the gaps mea- 
sured in thin films of Bi2Se3 on SiC substrate [23]. The 
experimental magnitudes are somewhat smaller which 
can be attributed to enhanced screening due to the pres- 
ence of the substrate which is neglected in our calcula- 
tions. The changes of band gaps are less systematic in the 
case of Bi2Te3. For slab thickness larger than 4QL, the 
magnitudes of Eg are larger in the calculations includ- 
ing the quasiparticle corrections. Second, the positions 
of Dirac points relative to the bulk bands change signifi- 
cantly. For both materials, the quasiparticle corrections 
"lift" the Dirac point from the bulk valence band. By 
extrapolating the results of our calculations on Bi2Se3 to 
the infinite slab thickness, we find that the incorpora- 
tion of quasiparticle corrections via the "scissors opera- 
tor" technique changes the position of Dirac point from 
0.04 eV below the bulk valence band [22| to 0.07 eV above 
it. The latter value is in better agreement with experi- 
mental results of Analytis et ai: 0.205 eV below the bulk 
conduction band or 0.095 eV above the bulk valence band 
assuming a 0.30 eV band gap |24j] . For Bi2Te3, the quasi- 
particle corrections change the Dirac point energy from 
-0.20 eV to -0.10 eV relative to the bulk VBM. The 
experimentally observed value is —0.13 eV [lj|. On the 
contrary, we find that some properties are not affected 
by the quasiparticle corrections. For instance, the de- 
gree of spin polarization of surface states investigated in 
Ref. [22| change very little. 

In conclusion, while the Kohn-Sham DFT band struc- 
tures are able to provide qualitative description of the 
topologically nontrivial electronic structure of Bi2Se3 and 
Bi2Te3, a quantitative agreement with a number of exper- 
imentally measured properties is achieved only after in- 
cluding the GW quasiparticle self-energy corrections. We 
further propose an energy-dependent "scissors operator" 
technique which allows the introduction of parametrized 
quasiparticle corrections into standard DFT calculations 
before SO interactions are included, thus greatly enhanc- 
ing their predictive power in describing systems based on 
the discussed topological insulators. 
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